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Abstract. We present the discovery of a binary companion to the pulsating sdB Feige 48. Using HST/STIS and 
archival FUSE spectra, we measure a period of 0.376±0.003 d and a velocity semi-amplitude of 28.0±0.2 kms~^. 
This implies that the companion star must either be of very low mass, or the orbit is at high inclination. Combining 
2MASS fluxes, the lack of a reflection effect, results from asteroseismology and a measurement of the rotation 
velocity of Feige 48, we show that the orbital inclination must be <11.4° and that the unseen companion is a 
white dwarf with mass >O.46M0. The implications of this discovery, and of binarity amongst sdB pulsators, is 
then discussed in the context of recent theoretical work on sdB formation. In particular we suggest that radial 
velocity studies focus on sdB pulsators with no known companion, and that asteroseismological studies of sdBs 
investigate a larger mass range than previously considered in order to test formation models. 
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1. Introduction 

I In recent years, interest in hot subdwarf B (sdB) stars 
' has increased dramatically. This is in large part due to 
I the discovery of rapid oscillations in some sdBs, but also 

■ because improving instrumentation has allowed more de- 

■ tailed studies of their spectral properties. It is now gen- 
' erally accepted that sdBs are associated with stars on 
, the extreme Horiz ontal Branch (EHB) fe.g.. lHebeil ll986': 
' ISafFer et al.l 1199^1 . This means they are low-mass stars 
. (~0.5Mq) that burn helium in their cores and have hy- 
' drogen envelopes too thin (<O.O2M0) to sustain nuclear 

burnin g. Stellar evolution calculations by iDorman et alJ 
l)l993|) predict that because of their thin envelopes, sdBs 
will evolve directly to the white dwarf cooling track, by- 
passing the Asymptotic Giant Branch. 

While the next stages of sdB evolution appear to 
be known, the question of sdB formation remains unan- 
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swered. There have been three basic scenarios pro- 
posed: binary interaction involving Roche lobe over- 
flow ijMengel et al.l Il976|) . single star evolution with 
strong mass loss near the tip of the Red Giant Branch 
(D'Cruz ct al. 199j), and the m erger of two helium white 
dwarfs l|lben fc Tutukovl 1X986^ . Significant progess ha s 
been made with the recent work of lHan et alJ l)2002ll200^ . 
who studied three possible binary formation channels: 
common envelope ejection, stable Roche lobe overflow and 
the merger of two helium white dwarfs. Using binary pop- 
ulation synthesis, they created 12 simulations and found 
one of them could satifactorily reproduce the observed 
characteristics of sdBs (Tgfj, log g, binary fraction, orbital 
period distributio n, etc) when c omp ared to radial velocity 
studie s bvlMaxted et alJ 1 I2OOI') and ' Morales- Rueda et alJ 
Napiwotzki et alJ lf2004ni fou nd a binary f ractio n 
of 39%, lower than predicted, while iLisker et all l)2004j) . 
in the largest homogenou s study so f ar, co mpared their 
data to the simulations of iHan et alJ l)2003|) and found a 
different best-fit model. 

The oscillations detected in some sdBs offer a way to 
test these formation models. Importantly, pulsations are 
seen in sdBs with short- and long-period binary compan- 
ions, as well as in apparently single sdBs. The properties 
of the V361 Hya stars (the short-period, p-mode sdB pul- 
sators, or i ginally called EC 14026 st a rs) ar e reviewed by 
iKilkennvl l|2002|i . while iGreen et all l)2003|) presents the 
detection of long-period, likely g-mode pulsations in some- 
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what cooler sdBs (currently known as the PG 1716+426 
or "Betsy" stars). 

Pulsat ions in the y361 Hy a star, Feige 48, were first de- 
tected bv lKoen et alJ l|l998t ). It is one of the coolest sdB 
pulsators in its cla ss, with Teff=29 500K and log 5=5.50 
llHeber et al1l200f]l) . and was discovered to be pulsating 
with five periods between 344 and 379 s. The oscillation 
amplitudes show dramatic variation from night to night, 
which the authors suggested are due to either changes 
in pulsation power, the beating of unresolved modes, or 
a combina tion of the two. In a comprehensive follow- 
up study, iReed et al.l l)2004j) observed Feig e 48 over 5 
years, and found that the periods measured bv lKoen et alJ 
l|1998^ were resolved and that the amplitudes vary by at 
least 30%. A tentative model match to four out of five de- 
tected periods, suggests that Feige 48 rotates with a period 
around 10 h. The Tgff and logg of this model were cho- 
sen to be consistent with those derive d from spectroscopy. 
Long-term monitoring bv lReed et alJ also provided an op- 
portunity to measure the phase stability of the pulsa- 
tions. In particular, they were able to place constraints 
on the period and mass of a companion to quite high 
precision, limiting any stellar-mass companion to have 
a period shorter than 3 days. In another recent study, 
ICharninet et alJ l(2003t) found a very good model match 
to all nine periods detected by high S/N photometry us- 
ing the CFHT. They de rived a rotat ional period of 9.58 h, 
in good agreement with lReed et al.l 

In this paper we present the detection of binary mo- 
tion in Feige 48, an sdB that was generally assumed to 
be single. We discuss the nature of the companion, and 
implications of this detection for current theories of sdB 
formation and pulsation. 



2. Observations 

Observations with HST/STIS were made as part of an 
ongoing project to measure metal abundances in pul- 
sators and non-pulsators (proposal #8635, PI: U. Heber). 
Analysis of these data will b e presented in a future pa- 
per fsee lO'Toole et al"l l2003[ for preliminary results). 
Once velocity variations were detected in these spectra, 
we searched both the HST and FUSE archives looking 
for any other observations. We found two sets of FUSE 
observations of Feige 48 in the archive, the first (proposal 
no. B033, PI: G. Fontaine) was a companion project to our 
HST/STIS observations, and their goal was also to mea- 
sure metal ab undances. The first re sults from this analysis 
is presented in lChaver et al. I ll2fl0,1) . The second, larger set 
of data was part of a project to study molecular gas in the 
interface between the galactic disk and halo (proposal no. 
C036, PI: R. Benjamin). Results from t his investigation 
will b e presented in a forthcoming paper ijBeniamin et alJ 
l2004l^ . Fortunately for us, the star used to illuminate the 
gas cloud under investigation was Feige 48. The list of 
observations are shown in Table ^ and their description 
follows below. 



Table 1. HST and FUSE observation log for Feige 48. 



Date 


UT 


exp. time (s) 


Instrument 


2001-03-26 


18:57:38 


2662.828 


FUSE 


2001-03-26 


20:37:27 


2567.906 


FUSE 


2001-05-13 


06:53:42 


3600.200 


HST/STIS 


2001-05-13 


07:59:32 


3600.200 


HST/STIS 


2001-05-13 


09:09:23 


3500.199 


HST/STIS 


2001-05-13 


10:13:33 


3511.199 


HST/STIS 


2002-04-03 


19:52:36 


3856.992 


FUSE 


2002-04-03 


21:29:57 


4226.070 


FUSE 


2002-04-03 


23:14:48 


4181.148 


FUSE 


2002-04-04 


00:59:19 


4264.227 


FUSE 


2002-04-04 


02:43:37 


4379.312 


FUSE 


total: 




41175.457 





2.1. HST/STIS UV echelle spectra 

We observed 5 sdBs with HST/STIS: three pulsators 
(Feige 48, PG 1219-^534 and PG 1605-^072) and two non 
pulsators (Feige 66 and Ton S-227) . The pulsators were all 
observed in time-tag mode, made possible by the use of the 
MAMA detectors. This was not originally intended, but 
was suggested (fortunately) by the time allocation com- 
mittee. Each star was observed using two medium reso- 
lution echelle gratings, E140M and E230M, giving useful 
wavelength ranges of 1150-1720 A and 1630-2370A respec- 
tively. The slit length and width used for all exposures 
was 0.2" x0.06", leading to spectral resolutions of 0.017 A 
at 1200 A and 0.020 A at 1500 A with the E140M grat- 
ing and 0.051 A at 1700 A and 0.068 A at 2400 A with the 
E230M grating. All spectra were processed using the time- 
tag and standard reduction routines in the STSDAS package 
in IRAF. The orders containing Lya were excluded, since 
they are not useful for a radial velocity analysis. 

2.2. FUSE spectra 

In both sets of observations, Feige 48 was observed using 
the high-throughput 30" x 30" LWRS aperture on all four 
spectrograph channels and recorded in time-tag mode. All 
data were reprocessed by the CalFUSE pipeline software 
(version 2.4). Of the four channels available, only the LiF2 
channel was used, since it contains no Lyman lines which 
are affected by interstellar absorption and airglow lines 
from the Earth. The detector also has the highest sensi- 
tivity of all onboard FUSE. The dispersion of this chan- 
nel is 6.7mA pixel^^ and the wavelength range covered is 
1086-1182 A. For a description of the detectors and mor e 
information on the FUSE mission, see IMoos et al.l()2000|) . 

3. Synthetic Spectra 

The UV spectra of Feige 48 (and sdBs in general) contain 
a large number of spectral lines due to iron-group elements 
(V, Cr, Mn, Fe, Co, Ni), making continuum definition of- 
ten difficult. The number of lines complicates abundance 
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Fig. 1. Velocity variations from far and near UV spectra 
of Feige 48. 

analyses, but is very advantageous for measuring veloci- 
ties, since they are very sharp (rotation velocities in sdBs 
are typically <5kms"^). 

As input for our spectrum synthesis, we used a metal 
line-blanketed LTE model atmosphere with solar metallic- 
ity and Kurucz' ATLAS6 Opacity Distribution Functions. 
The temperature, surface gravity, and helium abundance 
were set at r,ff=29 50 0K, log:q=5 50, an d log(He/H)=- 
2.93, as determined bv lHeber et"al] l|200(i> . 

The spectra were synthesised using Michael Lemke's 
version of the LINFOR program (originally developed 
by Holweger, Steffen, and Steenbock at Kiel University). 
Oscillator strengths were taken from the Kurucz line list, 
as were damping constants for all metal lines. Both the ro- 
tation and microturbulent velocities were set to Okms^^. 
(The rotation velocity of Feige 48 will be discussed further 
in Section O. 

4. Binary motion 

During the analysis of abundances in Feige 48, we found 
that the spectra from the two NUV orbits showed different 
velocity shifts compared to our synthetic spectra. Because 
the spectra were taken in time-tag mode, we were able to 
divide each of the far and near UV exposures into three 
1150 s sub-exposures. These spectra were cross correlated 
with our synthetic spectra and the results are shown in 
Figure^] The original FUSE observations were also taken 
in time-tag mode, allowing a similar division, with the 
exposure time typically being 1050 s. 

Using the double-precision cross-correlation routine 
(rvx) available in iraf,^ we measured the Doppler shift 
of each newly-created spectrum relative to our synthetic 
spectra. Since the HST/STIS spectra are in echelle for- 
mat, each observation yielded 40 velocity measurements 
(one per order). We calculated the weighted average of 
these velocities using the inverse square of the velocity 

^ available from http : / /iraf . noao . edu/extern/ 



Fig. 2. Amplitude spectrum of the radial velocities of 
Feige 48. The inset shows the aliasing around the adopted 
period. 

errors as weights. The end result of this was a list of ve- 
locities with two large gaps in time. These velocities are 
typically accurate to ~lkms^^. 

The analysis of the radial velocity curve is made sim- 
pler because the system is single-lined. By fitting the curve 
with a range of periods we produced a "power spectrum" , 
shown in Figure |21 which shows the goodness of fit. The 
actual period described below was found by selecting the 
highest peak in this spectrum. There is strong aliasing, 
leading to a relatively large period error (derived by mea- 
suring the difference between the nearest alias peaks). 

The ephemeris of the system at time To, defined as the 
conjunction time that the sdB moves from the blue to red 
side of the velocity curve, is 

HJD(To) = (2 452 367.05 ±0.03) 

+ (0.376 ± 0.0003 ± 0.003) x E 

where the first period error quoted is derived from the un- 
certainty due to the timespan of the observations, while 
the second error is derived by measuring the difference 
between the nearest alias peaks, as discussed above. We 
stress here that these values should not be taken as la 
uncertainties. We measure the velocity semi-amplitude 
of Feige 48 to be i4:sdB=28.0±0.2kms-\ while the sys- 
tem velocity is 7=-47.9±0.1kms^^. The velocity curve 
phased with this ephemeris is shown in Figure |2I Despite 
the separation in time of the observations, the phase 
coverage is quite good. The period and semi-amplitude 
of the sdB can be used to derive the mass function, 
/(m)=0.000856±0.000019MQ. If we assume a mass of 
^0.47 Mq for the sdB, as estimated from asteroseismol- 
ogy, then the lower limit for the companion mass is 
0.062 Af0. Note that the period we have derived here is 
well below the up per limit set using the phase stability of 
the pulsations bv lReed et al.l l)2004|) . Their analysis is in- 
sensitive to systems with short periods, small inclinations 
and massive companions, however, since the light-travel 
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Fig. 3. The phased radial velocity curve of Feige 48. The errors shown are formal errors determined by iraf. 
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Fig. 4. Spectral energy distribution of Feige 48. The filled 
diamonds are Johnson magnitudes, including 2MASS 
magnitudes. The HST/STIS and FUSE spectra are also 
shown. There is no evidence for the companion to the sdB. 

time across the system is far too small to cause a mea- 
sureable effect on the pulsations. 

5. The nature of the companion 

In order to discuss the i mplications of th is discovery for 
the formation models of i Han et alJ l)2003l) . we must first 
constrain the nature of the companion to Feige 48. 

The small mass function implies three possible com- 
panion types for Feige 48: a late-type main-sequence star, 
a brown dwarf or a white dwarf. We can place lim- 
its on the first two types using multi colour photome- 
try, as was done by iHeber et alJ l|2003 ^ for the sdB bi- 
nary HD 188112. Previous studies have suggested that any 



main sequenc e companion mu st have spectral type later 
than mid K (:,Koen et alJIlQQsh . We can now compare the 
2MASS JHK magnitudes, along with other Johnson mea- 
surements available, with the model spectral energy distri- 
bution of Feige 48, shown with photometric measurements 
overplotted in Figure^) Apart from 2MASS, the photom- 
etry is obtained from a variet y of sources (U — B = —1.03 . 
B -V = -0.25, V = 13.48 - iBern & WramdemarklHoH 
R = 13.57 (averaged), I = 13.72 - the USNO-B cata- 
logue, [Mone^^D'ioOS^. Magnitudes were converted to 
fluxes using th e calibration of Fukugita et al.. 1,1995,) for 
the optical and lCohen et al.| l)2003() for the infrared. The 
model was then scaled to the V flux of Feige 48, implying 
a distance to the system of 755±90pc. As can be seen in 
Figure^ there is no contribution to the flux from the com- 
panion in the if-band (on this scale, the errors in JHK 
magnitudes are smaller than the data points). Using the 
absolute K ma gnitude of dwarfs of spec tral types M2-M9 
determined bv^Beuer mann et al.l l|l999() . we find that if 
the companion to Feige 48 is on the main sequence, its 
spectral type must be M3 or later. This type of star will 
have a mass of ~0.35Mq, leading to a lower limit for the 
inclination angle of ^^14°. 

We can also consider information that has b e en ob - 
tained from the pulsati ons themselve s. 'Reed' 
iReed et all ()2004j) and .Charoinet et all (i2003) have car- 
ried out asteroseismological analyses of Feige 48, and sug- 
gest a rotation period of ^9.5-10 h. Given the uncertainties 
involved in the sdB evolution models used and the aliasing 
seen in Figure El this is in very good agreement with the 
orbital period of 9.03 h determined in Scction|31 In partic- 
ular, both .Reed [ 2001) and Reed et al. ( 2004 ) assumed a 
rotational splitting coefficient (C„i) of zero. If C„; is set to 
0.1, which is not unreasonable for a slightly evolved sdB, 
then the rotation period would be ~9 h. 
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1.2 




Fig. 5. Fit of V sin i to three Fe III Unes. The Une at ^^2098. 5 A is blended with a Mn ill hne. The dashed line represents 
wsini=Okms~"'^, the dotted line is wsini=5kms^^, and the dash-triplc-dot line is wsin j=10kms~^. See the text for 
further details. 



We now compare the observed upper limit on a re- 
flection effect to what might be expected from an M3 
dwarf with radius ~O.35i?0. A similar system, but with 
higher inclination and shorter period is the sdB+dM 
HS 2333+3927 ( Heb er et al 2004). In this case a large 
reflection effect is observed (Ay=0.28mag) even though 
the inclination is only 36-39°. This suggest that an effect 
might be measureable for Feige 48 if present. Calculations 
by H. Drechsel (private communication) find that an 
sdB+dM system with a period of 9 h and inclination 15° 
would show luminosity variations of 0.025-0.030 mag for 
maximal albedo of the secondary. The only time-series ob- 
servations of Feige 48 made were optimised for oscillation 
studies, not searching for a reflection effect. With an or- 
bital period of 9 h, such an effect is difficult to detect, since 
observations from a single site will typically not cover an 
entire orbit. A re-analysis of these data places an upper 
limit of '^O. 013 mag on any reflection effect (M. Reed, pri- 
vate communication). Given the uncertainties involved in 
these calculations and comparisons, we do not feel con- 
fident completely ruling out an M3 dwarf companion at 
this stage, but instead suggest the white dwarf scenario is 
more likely. 

Finally, using our high resolution UV spectra, we may 
be able to improve on the equ atorial rotation velocity lim- 
its set bv lHeber et all l|200n|) and therefore constrain the 



orbital inclination and companion mass. In the near UV 
our spectral resolution i s almost a fa ctor of two better 
than the spectra used bv lHeber et all so we shifted each 
of the six spectra used to determine radial velocities to 
their rest wavelength and then averaged them. The spec- 
trum was examined for metal lines in close blends that are 
clearly resolved (a difficult task due to the very large num- 
ber of unresolved blends), and these lines were fit with a 
variety of values of wsini. We show the three Feiii lines 
chosen in Figure [3 with three different rotation velocities 
wsini=0, 5, and lOkms"^. There is no clear difference 
between the Okms^^ (dashed) and 5kms~^ curves (dot- 
ted), while wsini = lOkms^^ (triple-dash-dotted curve) 
does not appear to fit. After a simple analysis, we find 
that wsini = 4.8 ± 1.2kms^^. At a 3a level usini is con- 
strained to less than 8kms~^. Together with the upper 
limit of 5kms~^ (3(t) der ived from the Mgii line in the 
optical ifHeher et a,lJl200^^ . we consider 5 km s as an up- 
per bound the following discussion. 

The rotation of Feige 48 is expected to be tidally locked 
to the orbit. In this case we can estimate the rotation ve- 
locity from the orbital period and the radius (^ 0.2i?Q) 
to be 26.9li[ g km (this is coincidcntally approximately 
the same as the orbital velocity). Comparison with the 
projected rotational velocity constrained above yields an 
upper limit to the inclination of the system of i = 11.4°. 
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Thi s is the same upper limit derived bv lReed et alJ l)2004|) 
and iReedl l)200ll) when they com bined the upper lim it to 
the rotation velocity derived by iHeber et al.l l)200Cl|) and 
the rotation period derived from their asteroseismologi- 
cal analysis. From the mass function we derive a lower 
limit to the companion mass of M2 = 0.46Mq. This re- 
sult means that the companion must be degenerate, and 
is most likely a white dwarf. Considering the orbital pe- 
riod and companion mass, the Feige 48 system should have 
gone through two common envelope phases, and therefore 
can be explained by the " second common e nvelope ejec- 
tion channel" described bv lHan et alJ l)2003j) . 

6. Discussion and Conclusions 

We have detected binary motion in the V361Hya star 
Feige 48 with a period of 0.376 d and a velocity semi- 
amplitude of 28.0kms^^. Combining multicolour photom- 
etry, the asteroseismological determination of the rotation 
period, and an upper limit of the equatorial rotation ve- 
locity, we conclude that the orbital inclination i <11.4°. 
This implies that the minimum mass of the companion is 
0.46 Mq; combining this with the lack of a reflection effect 
suggests the star must be a white dwarf. 

Feige 48 is one of only two pulsating sdBs with a 
known white dwarf companion. Unlike KPD 1930-1-2752, 
the other system, it has a relatively simple pulsation 
spectrum and no (or v e ry lit tl e) ellipsoida l disto rtion. 
This has allowed iReedl l|200l|) . iReed et al.l (|2004|) and 
ICharpinet et alJ l(2003^ to match pulsation frequencies 
and constrain sdB evolution models. However, these stud- 
ies assumed that Feige 48 evolved according to canonical 
EHB theory. While this is probably accurate since the 
star has undergone two p hases of co mmon envelope evo- 
lution, the conclusions of iHan et alJ (|2003) suggest that 
investigating a wider mass range than previously may 
be warranted, with masses as low as ~O.36M0. In the 
case of single sdBs, which should be formed through the 
merger channel, masses in the range 0.4-0.7 Mq should be 
examined. The sdBs in long-period orbits with main se- 
quence companions (formed through the Roche lobe over- 
flow channel) may fall in t he 0.3-0.5 M m range, and pos- 
sibly even higher, although iHan et aD find this to be un- 
likely. Enlarging the mass range investigated in this way 
may provide another useful test of the proposed formation 
models for sdBs. 

These tests can only be carried out once the binary 
nature of the pulsating sdBs is known. Up to now only 
seven V361 Hya stars are confirmed binaries; four have 
F- or G-type main sequence companions (V361 Hya itself, 
PB8783, EC 20117-4014 and EC 10228-0905), two have 
white dwarf companions (KPD 1930-1-2752 and Feige 48), 
and one is in an eclipsing binary system with an M5 main 
sequence star (PC 1336-018). That leaves around 25 pul- 
sators with unknown status. The fraction of sdBs is binary 
systems is still unde r debate, howeve r if we consider the 
best-fit simulation ofl Han et alJl|2003l) . we find that '^41% 
or 10-11 of these stars should have either late-type main 



sequence or white dwarf companions. iLisker et all l^2004^ 
compared their observations of 52 sdB stars wit h the 12 
simulation sets constructed by Han et aD l|2003l) using a 
simple statistical test. The simulation set they found to 
fit their obser vations differs from that found qualitatively 
bv lHan et al.l and suggests that up to ~53% or ~13 of the 
pulsating sdBs should have an unseen companion. High 
precision spectra are needed to observe short period sys- 
tems such as Feige 4^8, so a survey with high spectral and 
temporal resolution is required. 
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